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Abstract 

We consider a supersymmetric model in which gaugino masses are generated by 
the anomaly-mediation mechanism while scalar masses are from tree-level supergravity 
interaction. In such a model, scalar fermions as well as Higgsinos become as heavy as 
0(10—100 TeV) and hence only the gauginos are superparticles kinematically accessible 
to the LHC. We study how and how well the properties of gauginos can be studied. 
We also discuss the strategy to test the anomaly- mediation model at the LHC. 



1 Introduction 



Anomaly mediation of supersymmetry (SUSY) breaking in a hidden to the SUSY standard- 
model (SUSY SM, or SSM) sector [H [2] is very attractive, since it is the simplest mechanism 
for the mediation of SUSY breaking. Namely, the anomaly mediation always takes place in 
a generic SUSY theory and hence we do not need any extra assumption beside the presence 
of SUSY breaking sector to mediate the SUSY-breaking effects to the SSM sector. Without 
any special requirement on the Kahler potential, but just by assuming that there is no singlet 
field in the SUSY breaking sector, the scenario of anomaly mediation predicts the so-called 
split SUSY spectrum [3] where squarks and sleptons may have masses of the order 10 — 100 
TeV while the masses of gauginos are in the range of 100 GeV — 1 TeV |4, 5jcl 

Because of the relatively large masses of squarks we need a very precious fine-tuning of 
parameters to obtain the correct electroweak symmetry breaking, but on the other hand it 
solves many serious problems in the SSM. First of all the flavor-changing neutral current and 
CP-violation problems become very milder due to the large masses of squarks and sleptons. 
We may naturally explain no discovery of Higgs at LEP and no discovery of proton decays 
induced by dimension-five operators. Furthermore, the gravitino mass is also predicted at 
the order of 100 TeV, which makes the cosmological gravitino problems much less severe 
[HI E]. In fact, it has been pointed out that the leptogenesis [8j does work in the anomaly- 
mediation model, since the reheating temperature Tr can be as high as 10^° GeV without 
any conflict with cosmology p]. 

Although the squarks and sleptons are so heavy, the masses of gauginos may be in the 
accessible range to the LHC experiments. Thus, even in this model, discovery of the signals 
from the productions of superparticles may be possible at the LHC. More importantly, the 
anomaly-mediation model predicts unique mass relation among gauginos, as we will see 
in the following. In particular, the anomaly mediation predicts m^r < nij^ < nig (with 
m^, m^y, and rrig being the gaugino masses of f/(l)y, SU(2)l, and SU{3)c gauge groups, 
respectively) in a large region of the parameter space [5]. In this letter, we discuss how the 
anomaly-mediation model can be studied at the LHC. 

2 Model 

Let us start our discussion by summarizing basic features of the model. As we have mentioned 
in the previous section, we consider the anomaly-mediation model in which all the sfermions 
as well as heavy Higgses and Higgsinos are heavy (of masses of the order of 100 TeV). 
Assuming that there is no singlet field in the SUSY breaking sector, all the gaugino masses 
are suppressed compared to the SUSY breaking scalar masses. 

In this class of models, gaugino masses are mainly from the effect of anomaly mediation 

T^^In this letter, we use "anomaly-mediation model" for those where gaugino masses are generated by the 
effect of anomaly mediation while the scalar masses are from tree-level supergravity interaction between 
observable-sector fields and SUSY breaking fields. 
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and possible radiative correction due to the Higgs-Higgsino loop [H [TO]. Then, the SUSY 
breaking gaugino mass parameters (at the scale of sfermion masses mj) are obtained as 



M, = jL(m3/2 + L), (2) 
Ms = j|^(-3m3/2), (3) 

where gi, g2, and are gauge coupling constants of U{1)y, SU{2)l, and SU{3)c gauge 
groups, respectively, and 7/13/2 is the gravitino mass. (Hereafter, we use the convention such 
that m3/2 is real and positive.) In addition, 

2 I 1 2 

L = fism2P—P^\n^^, (4) 
l/ip — 

with tan jS being the ratio of two Higgs bosons, and the mass of heavy Higgses. Thus, 
due to the Higgs-Higgsino loop contribution, gaugino masses may significantly deviate from 
the pure-anomaly- mediation relation (which is given by taking L = 0). However, since the 
natural sizes of and ttla are both of the order of the gravitino mass, the L-parameter 
is expected to be at most 0(ms/2)- (Notice that \L\ < m^.) Consequently, with a natural 
choice of fi- and 5^-parameters, Wino becomes the lightest among gauginos in a large region 
of the parameter space. In this case, taking account of the radiative correction due to 
electroweak gauge bosons, neutral Wino becomes the lightest superparticle (LSP) [llj. 
Charged Wino becomes slightly heavier than the neutral one; the mass splitting is as 
small as — m^^o ~ 155 — 170 MeV. As we will see, the smallness of the mass splitting 
has an important implication to the study of the anomaly-mediation model at the LHC. 

We should also note that three gaugino masses depend on three parameters: ms/2, \L\, 
and Arg(L). Even so, non-trivial constraint exists among the gaugino masses. In order to 
see how the gaugino masses are constrained, it is instructive to approximate the physical 
(i.e., on-shell) gaugino masses by Mi — Ms given in Eqs. ([I]) — ([3]). Then, we can see 
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Thus, once the Wino and gluino masses are fixed, upper and lower bounds on the Bino mass 
are obtained. (In our quantitative analysis in the following sections, we calculate the bounds 
more accurately by taking into account the renormalization-group effect below the sfermion- 
mass scale m^, taking ttl^ = We have checked that the dependence on is rather 

mild; the bounds change ~ 10 GeV when is varied in the range ms/2/2 < rrij < 2m^i2-) 
It is an important test of the anomaly-mediation model to see if the gauginos satisfy the 
mass relation. 
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Figure 1: Cross section for the process pp — gg for ^/s = 14 TeV. Here the squarks are 
assumed to be heavy enough to be neglected. 

3 Anomaly- Mediation Model at the LHC 
3.1 Set up 

Because the gauginos are the only superparticles kinematically accessible to the LHC, and 
also because an important information is imprinted in the gaugino masses, we must study the 
properties of the gauginos at the LHC if the anomaly-mediation model is realized in nature. 
In the following, we discuss how well the gauginos can be investigated at the LHC. As an 
example, in this letter, we consider the case where the underlying model gives fn^ = 400 GeV, 
m^y = 200 GeV, and rrig = 1 TeV (which is given by 171^/2 — 39 TeV, \L\ ~ 28 TeV, and 
Arg(L) = 0). Notice that the neutral Wino is the LSP. 

First, we comment on the decay of the charged Wino. Since the mass splitting m^y± —m^^o 
is very small, decays into W'^ emitting an extremely soft pion. Thus, it is not easy to 
identify the decay of the charged Wino at the LHC. In addition, since the typical decay length 
of is 0(1 — 10 cm), it is highly challenging to find the track of the charged Wino although 
it may not be impossible. Thus, in our main discussion, we make a conservative assumption 
that we cannot find the charged- Wino track; then both and are treated as invisible 
particles. However, we will also discuss the implications of discovering the charged- Wino 
tracks. 

The most important production process of superparticles at the LHC is the pair produc- 
tion of gluinos ^'s: pp gg, as far as the gluino mass is about 1 TeV or smaller; in Fig. 
[1] we plot the cross section for this process. When rrig ~ 1 TeV, for example, thousands 
of supersymmetric events will be available at the LHC with C = 10 fb"^ (with C being 
luminosity) . 

Once gluino is produced, it decays into a lighter gaugino (i.e., Wino W or Bino B) 
and standard-model fermions. If Bino B is produced by the decay of gluino, it decays 
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successively to Wino. Experimental signals strongly depend on how the gauginos decay. 
The decay patterns of gauginos are sensitive to the masses of sfermions and Higgsinos, and 
hence are not predictable. Thus, we make several assumptions for our study. For the decay 
of gluino, we adopt 

Br{g Bqq) = 1 - Br{g Wqq) = 0.25, (6) 

where q denotes the standard-model quarks. 

On the decay of Bino, important processes are B W^W^ via Bino-Wino mixing, 
B — > W^hsM (with hsM being standard-model-like Higgs boson) via Bino-Higgsino mixing, 
and B Wff (with / being standard-model fermions) due to sfermion-exchanges. In a 
large fraction of the parameter space, B W^W"^ and B W^h^u become the dominant 



decay modes|^ If so, however, it is difficult to study the properties of B. Then, probably 
the anomaly-mediation model can be tested only by finding the charged- Wino track, which 
will be challenging, (as well as by the measurement of rrig — m^y, which we discuss below). 
If there exists a significant hierarchy between sfermion and Higgsino masses, however, the 
last process may dominate. In particular, if the sleptons are much lighter than Higgsino 
(by factor 10 or so), B —>■ Wl^l^ may acquire a large branching ratio. As we will see, 
invariant-mass distribution of gives an important information on the Bino mass. Thus, 
in order to demonstrate what we can learn from the decay mode B — >■ W^l^l~ , we consider 
a special case where sleptons are much lighter than the Higgsino so that the decay modes, 
B — i> W'^W'^ and B W^hsu, become negligible. For our simulation, we use 

Br{i3 WLL) = 1 - 5r(5 ^ Wqq) = 0.3, (7) 

where L = is for the standard-model leptons. 



3.2 Simulated background samples and selections 

The gluino pair is produced mainly with gluon-gluon collision, and the gluino decays into 
W oi B with two or four high pt jets. Therefore the missing transverse energy, , carried 
away by two Winos plus four high pt jets is the leading experimental signature (defined as 
"no lepton mode"). Also lepton pair from the decays of B, four-jet and^^ is the next leading 
signature of the signal (defined as "dilepton mode"). Fraction of the event with dilepton 
depends strongly on the decay branching fraction, Br{cj Bqq) x Br{B W^l~^l~). The 
discovery potential will be mainly determined with the no lepton mode and the dilepton 
mode will be used to measure the mass difference between W and B. 

The following four standard-model processes can potentially have the event topology of 
I^x with jets: 

• + jets, {W^ (u) 

• + jets, (ZO uu,T+T-) 

^^We thank T. Watari for a useful discussion on the Bino decay. 
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Figure 2: distribution of the SUSY signal and background processes: The open histogram 
shows the SUSY signal (with = 1 TeV) and the hatched shows the sum of background 
distributions. 

• tt 

• QCD jets: Heavy flavor quarks (6 and c) with semi-leptonic decay and the light flavor 
jets with mis- measurement. 

The high pt multi-jets are key of the analysis and they should be estimated with the 
Matrix-Element calculation. Parton-Shower is not good approximation in such a high re- 
gion and the background contributions would be underestimated [12j . These background pro- 
cesses are generated with ALPGEN2.05 [13], and the exact Matrix- Element Calculations are 
applied up-to five partons. The produced 200 M events are fed into the Parton-Shower gen- 
erator (JIMMY4.0/Herwig6.5 [H]) in order to evolute the QCD shower. Multiple-interaction 
processes are also taken into account. Matching between Matrix- Element and the evaluated 
Parton-Shower are also applied to remove the double counts. The detector effect is taken into 
account using the smearing Monte Carlo simulation of the ATLAS detector (ATLFAST [15j). 

The following simple event selections are applied, 

• Number of jets with pt> 50 GeV is larger than or equal to 4. 

• pt of the jets are required to be larger than 200 and 100 GeV for the leading and 2nd 
leading jets, respectively. 

• The missing transverse energy, , is larger than 300 GeV. Fig. [2] shows the 
distributions for the no lepton mode. 

• The same flavor two leptons with pt> 20 GeV are required for the dilepton mode. 

The effective mass, which is define as I^t +J24jets Pt, is a good variable to discriminate 
the SUSY signal from the SM background processes, and Figs. [3] show the effective mass 
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Figure 3: Effective mass distributions of tlie SUSY signal and background processes: No 
lepton mode and dilepton mode. In both figures, tlie open liistogram shows the SUSY signal 
and the hatched shows the sum of background distributions. (Blue circle, red triangle, green 
triangle and magenta box show the top, and QCD processes, respectively.) 

distributions of the signal mentioned above and the SM background processes for the no 
lepton and dilepton modes. Three standard model processes, + jets, + jets, and ti 
contribute equally to no lepton mode, but the signal excesses (with rrig = 1 TeV) will be 
observed in the high effective mass region. ATLAS has a discovery potential up-to 1.2 TeV 
with an integrated luminosity of 10 fb~^. ti process is the dominant background for two 
lepton mode. Since flavor of the selected two leptons are independent in this case, these 
background processes can be removed easily using the flavor subtraction [TB] . 

3.3 Mass reconstruction 

Now, let us discuss how and how well the gaugino masses can be reconstructed. With the 
gluino pair production, we obtain, at the parton level, (i) four quark jets, and (ii) several 
charged leptons (or tau jets). The momentum distribution of the quark jets contains an 
information on and rrig, while that of charged leptons contains that on and m^. 
Indeed, if we concentrate on the decay of single gluino: g —>■ Wqq, invariant mass of the qq 
system Mgg is required to be 

Mgg <mg- Ul^ , (S) 

while, for the case of ^ — Bqq, followed by -B ^ W^l^l^, we obtain 

Mi+i- <m^- m^^, Mi+i-gg < m-g - ruy^,, (9) 

where M;+;- and Mi+i-^g are invariant masses of qq and l^l~qq systems, respectively. Thus, 
from the invariant- mass distributions, we can extract the information on the gaugino masses. 

The same selections mentioned in the previous subsection are applied to make the mass 
distributions for M^g, Mi+i- and Mi+i-qq. The integrated luminosity of 100 fb^^ is assumed 
to measure these invariant masses. The leading four jets are used to make combination of 
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Mjj [GeV] 

Figure 4: Mjj distribution for the no lepton mode. 



dijets, and there are three possible combinations of four jets, (1+3, 2+4), (1+4, 2+3) and 
(1+2, 3+4). For example, (1+3, 2+4) means that the leading and 3rd leading jets are 
considered to be emitted form the same gluino, and 2nd and 4th leading pt jets have the 
same parent. The fractions of the correct combination are the same in the combinations 
of (1+3, 2+4) and (1+4, 2+3). We select the combination in which the difference between 
two calculated invariant masses of dijets is smaller than the other combination. On the 
other hand, the fraction of the wrong combination is larger in the (1+2, 3+4) combination, 
so this combination is used for comparison with the other combinations. The difference 
between two calculated invariant masses is required to be smaller than 100 GeV in oder 
to suppress the contamination of g Bqq — > qqqqW . Fig. H] shows the distribution of the 
invariant mass of dijets Mjj, in which the fitted endpoint using linear function is 784+37 GeV 
(the expected endpoint is 800 GeV). Statistical error of the fitted endpoint is 5% with an 
integrated luminosity of 100 fb~^. Jet is reconstructed with the cone algorithm (size of 0.4), 
and some part of the evaluated shower escapes from the cone. This is the reason why the 
reconstructed Mjj becomes smaller than parton level (M^g). This shift can be collected 



with parton-jet calibrations and the uncertainty in this calibration is finally 1%|^ There 
is also another systematic uncertainty in the fitting procedure. The other samples with the 
different mass combinations are also generated, and the same selections and fitting procedure 
are applied. The fitted endpoints are proportional to the mass difference at parton level and 
consistent within an error of 5%. 

Fig. El^a) shows the Mi+i- distribution for the dilepton mode after the flavor subtrac- 
tion [T6] is performed. Clear edge is observed at the expected endpoint (200 GeV), and 
the statistical error of the fitted edge is 1 % with an integrated luminosity of 100 fb~^. 
Systematic uncertainty in the absolute energy calibration is much less than the statistical 



T^^The calibration of jet energy to parton level will be performed using W ^ qq and gq — *■ 75 (with g and 
7 being gluon and photon, respectively) events, and finally 1% level accuracy will be obtained at ATLAS 
detector. 
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Figure 5: Mi+i- and M;+;-qq distributions for the dilepton mode. 



error. 

In order to calculate the invariant mass of dijets and dileptons the corresponding two jets 
are selected as follows: One of the leading or 2nd leading jet is selected, which is closer 
to the dilepton system. The distance between the dilepton and this jet is also required to 
be smaller than 2.0. The 4th leading jet is considered as another jet, since the pt becomes 
smaller due to the smaller mass difference of g and B. The invariant mass between dilepton 
and the dijets are shown in Fig. [5](b), after the flavor subtraction is applied. Clear edge is 
still observed at 800 GeV, and the statistical error is 50 GeV. 

3.4 Other possible information 

So far, we have seen that some of the mass differences can be determined by the LHC. In 
particular, if Br{g Bqq) and Br{B —>■ W^l~^l~) are both sizable, two mass differences 
can be reconstructed. Even so, the confirmation of the anomaly-mediation mass spectrum 
of gauginos is not straightforward because the following questions, which are crucial to test 
the anomaly-mediation model, are not answered yet. 

First, we have seen that the information on mass differences is obtained from the invariant- 
mass distribution of jets and leptons. However, it is still an open question how the masses of 
gauginos can be determined. In addition, one should note that the event shape from the de- 
cay chain g — * Bqq followed by i? — ^ W'^l'^l^ (or B —>■ W^Pu) is hardly distinguished from 
that from g — > W^qq followed by 1^° Bl^l~ (or g W^qq' followed by Bl^u), 
in particular if the track of the charged Wino is not observable. Thus, it may not be easy 
to confirm that Wino is the LSP. For the test of anomaly-mediation model, it is important 
to answer the following questions: 

1. How do we determine the masses of gauginos rather than mass differences? 
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2. How do we know that the hghtest state is Wino, not Bino? 

These can be easily answered if the (short) tracks of charged Wino can be identified at 
the LHC [5]. Although it is very challenging, such tracks may be seen in the inner detector, 
like the semiconductor pixel detector and/or the transition radiation tracker (TRT). If the 
tracks of charged Winos can be identified, it is obvious that the lightest gaugino is Wino, 
not Bino. In addition, Wino mass may be determined by using timing information combined 
with momentum information on the charged- Wino track. The resolution of the determined 
P is about 0.1 if /3 is less than 0.85. So the mass can be determined with accuracy of 10%, 
if we have enough samples of this exotic track. 

Even if the information on the charged- Wino track is not available, the above questions 
may still be answered. For the first question, one way to determine the gaugino masses 
is to use the cross-section information on the process pp gg. As one can see in Fig. 
m a{pp gg) strongly depends on the gluino mass. Thus, even if the error in the cross 
section is fairly large, we may still have a useful constraint on the gluino mass. When the 
signal number can be determined very roughly with an accuracy of 20 % level, which is the 
uncertainty of the background, gluino mass can be estimated with an accuracy of 3 %. 

In order to discriminate the possibility of Bino-LSP, we may be able to use the Drell-Yan 
process pp — ^ WW. If the W^ is the lightest neutralino, such a process does not provide 
any signal. If Bino is the LSP, on the contrary, Winos produced by this process decay, 
resulting in events with large missing Et and possibly multi-leptons. From the negative 
search for such event, we may be able to discriminate the possibility of Bino-LSP. Another 
possibility is to tag the flavor of quarks emitted by the decay of gluino, which may be 
possible for third-generation quarks. In the Wino-LSP case, the following process may 
occur: g Bhh W^l^vhh while, for the Bino-LSP case, g W~^bi Bl^vhi. Thus, by 
identifying the first decay chain using 6-jet tagging, we may be able to confirm that the LSP 
is Wino. 

3.5 Testing the anomaly- mediation model 

Now we consider how well the anomaly-mediation model can be tested. With the measure- 
ments of gaugino masses discussed in the previous subsections, we may be able to check if the 
gaugino masses satisfy the mass relation in anomaly-mediation model discussed in Section 

El 

For this purpose, we first summarize the experimental constraints. For a given set of 
gluino masses, we plot the regions allowed by the expected LHC constraints on the 
vs. plane. Here, based on the discussion given in the previous section, we assume that 
— rriy^ and rrig — rriy^ can be experimentally determined with the errors of 1 % (i.e., 2 
GeV) from Mi+i- and 5 % (i.e., 40 GeV) from Mi+i-jj and Mjj, respectively. The result is 
shown in Fig. O (Notice that — m^, is determined so accurately that the upper and 
lower bounds on shown in the figure look almost like a single line.) On the same figure, 
we also show parameter region which is allowed in the anomaly-mediation model. (See the 
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Figure 6: Expected experimental constraints on m^^ vs. plane. The gluino mass is fixed 
to be 900 GeV, 1 TeV, and 1.1 TeV, from below. Solid lines (which are almost degenerate 
in the figure) are expected experimental upper and lower bounds on as functions of rriyy, 
while the dashed lines are upper and lower bounds on myy. In the shaded region, gaugino 
masses are consistent with the prediction of the anomaly-mediation model. The square in 
the middle figure shows the underlying gaugino masses used in our numerical analysis. In 
addition, with the gluino mass being fixed, the Wino and Bino masses can be calculated in 
the case of pure anomaly mediation; such point is shown by the star. 
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discussion at the end of Section [2l) One of the crucial test of the anomaly mediation is to see 
if those three regions overwrap for a gluino mass consistent with experimental constraints. 

The gluino mass should be constrained by using the cross section for the process pp — *■ 
gg, as discussed in the previous section. If we can experimentally confirm that the gluino 
production cross section is consistent with the predicted value for the case of rrig ~ 1 TeV, 
we can conclude that the gaugino mass spectrum is consistent with the prediction of the 
anomaly-mediation model. 

We also note that, in the present case, it is possible to exclude the pure anomaly- 
mediation model (i.e., L = 0). In the case of pure anomaly mediation, the Wino and 
Bino masses are determined once the gluino mass is fixed. Then, with the precise informa- 
tion on — rriyy from Mi+i-, prediction of the pure anomaly mediation can be excluded 
for the underlying parameter used in the current analysis. In general, accurate determi- 
nation of — will be very useful to exclude (or confirm) the prediction of the pure 
anomaly-mediation model. 

4 Summary 

In this letter, we have discussed how and how well the anomaly-mediation model can be 
studied at the LHC As we have seen, the masses of gluino and Wino may be constrained 
by using invariant-mass distribution of dijet as well as the cross-section information. 

On the contrary, the Bino mass is hardly studied in a large fraction of parameter space 
since the dominant decay modes of Bino are expected to be 5 — >• W'^W^ and B W^h^M- 
Therefore, it is very difficult to test the gaugino mass relation in the anomaly-mediation 
model at the LHC. However, in a special case where the sleptons are much lighter than the 
Higgsino, the decay mode B —>■ Wl^l~ may acquire a large branching ratio. In such a case, 
the mass difference between Bino and Wino can be well determined by using invariant-mass 
distribution of dilepton and hence the gaugino-mass relation may be tested. 

In any case, in order to confirm that Wino is the LSP, it will be important to find the 
charged- Wino track. Thus, if the excess of event is observed without the discovery of 
any sfermions, the search for the short-lived charged- Wino track (with the decay length of 
0(10 cm)) is strongly suggested. 

Note added in proof: The present procedure can be also used in models other than 
anomaly-mediation model, as far as all the sfermions are heavier than the gauginos. For 
example, one may test the grand-unified-theory relation among the gaugino masses. 

Acknowledgement: The authors would like to thank M. Ibe for useful discussion at the 
early stage of the project. We have used MadGraph/MadEvent packages for some of our 
numerical calculations. 

References 

[1] G. F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 9812 (1998) 027. 



11 



[2] L. Randall and R. Sundrum, Nucl. Phys. B 557 (1999) 79. 

[3] N. Arkani-Hamed and S. Dimopoulos, JHEP 0506 (2005) 073; G. F. Giudice and A. Ro- 
manino, Nucl. Phys. B 699 (2004) 65 [Erratum-ibid. B 706 (2005) 65]; N. Arkani- 
Hamed, S. Dimopoulos, G. F. Giudice and A. Romanino, Nucl. Phys. B 709 (2005) 
3. 

[4] J. D. Wells, Phys. Rev. D 71 (2005) 015013. 

[5] M. Ibe, T. Moroi and T. T. Yanagida, Phys. Lett. B 644 (2007) 355. 

[6] See, for example, M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 625 (2005) 7; 
Phys. Rev. D 71 (2005) 083502; K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 
73 (2006) 123511. 

[7] M. Kawasaki, F. Takahashi and T. T. Yanagida, Phys. Lett. B 638 (2006) 8; Phys. Rev. 
D 74 (2006) 043519. 

[8] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45; For a recent review, 
W. BuchmuUer, R. D. Peccei and T. Yanagida, Ann. Rev. Nucl. Part. Sci. 55 (2005) 
311. 

[9] M. Ibe, R. Kitano, H. Murayama and T. Yanagida, Phys. Rev. D 70 (2004) 075012. 

[10] T. Gherghetta, G. F. Giudice and J. D. Wells, Nucl. Phys. B 559 (1999) 27. 

[11] J. L. Feng, T. Moroi, L. Randall, M. Strassler and S. f. Su, Phys. Rev. Lett. 83 (1999) 
1731. 



[12] S. Asai and T. Sasaki, session 2.4 in hep-ph/0608322| 



[13] M.L. Mangano, M. Moretti, F. Piccinini, R. Pittau, A. Polosa, JHEP 0307 (2003) 001. 



14] G. Corcella et al., hep- ph/001 1363 http://projects.hepforge.org/jimmy/ 



[15] ATLAS Physics TDR vol.1 , CERN/LHCC/99-15. 

[16] ATLAS Physics TDR vol.2 , CERN/LHCC/99-15, p.825. 

[17] T. Stelzer and W. F. Long, Comput. Phys. Commun. 81 (1994) 357; F. Mahoni and 
T. Stelzer, JHEP 0302 (2003) 027. 



12 



